Boltzmann equation approach to transport in finite modular quantum systems 
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We investigate the transport behavior of finite modular quantum systems. Such systems have 
recently been analyzed by different methods. These approaches indicate diffusive behavior even and 
especially for finite systems. Inspired by these results we analyze analytically and numerically if 
and in which sense the dynamics of those systems are in agreement with an appropriate Boltzmann 
equation. We find that the transport behavior of a certain type of finite modular quantum systems 
may indeed be described in terms of a Boltzmann equation. However, the applicability of the 
Boltzmann equation appears to be rather limited to a very specific type of model. 
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I. INTRODUCTION 

There are essentially two major tools which are used 
to analyze the transport behavior of quantum systems: 
Linear response theory as implemented in terms of the 
Kubo formula 0, 0, & 0> 0> H > an d approaches based 
on the Boltzmann equation (BE) 0] . The validity of the 
latter has been subject to ongoing discussions during the 
last century Q . This refers to the BE as a method to de- 
scribe gas-dynamics on purely classical grounds. In the 
context of quantum mechanics the situation may even be 
more complicated. Can the dynamics of systems that are 
controlled by the Schrodinger equation (SE) be mapped 
on a BE? And if so, how? Considerable work in that 
direction has been done by Peierls @, Kadanoff, Baym 
[lfj| and others [TJ E2, EH- I n this literature it is fre- 
quently pointed out that the mapping typically relies on 
additional assumptions like the "random phase approxi- 
mation" [9( or the possibility to truncate graphic expan- 
sions that does not necessarily follow from the underlying 
dynamics [Iol |. But also recent publications address the 
mapping of quantum dynamics onto BE's [HI, EM EB] • In 
the article at hand we investigate the applicability of a 
BE to systems which are both, complex enough to exhibit 
diffusive behavior and simple enough to be analyzed from 
first principles by direct numerical integration. Thus the 
dynamics as resulting from a BE may simply be com- 
pared to the dynamics as resulting from the SE. The re- 
sults for transport behavior of those systems obtained by 
other methods (Kubo formalism, Hilbert space average 
method (HAM), time-convolutionless (TCL) [13, El, El, 
which have been mentioned in the abstract, may be found 
in @,[H,[23,[13 

The article at hand is organized as follows: 
In Sect. |TT] we very briefly review the concepts under- 
laying the famous BE which is meant to be a gross de- 
scription of the dynamics of dilute gases. We comment 
on linear forms of the BE and their diffusive solutions, 
i.e., we state an explicit form for the diffusion coefficient. 
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In Sect. IIIII we introduce our finite modular quantum 
system which may be viewed as model for a particle 
hopping between a few lattice sites or a model for the 
energydynamics within a chain of, e.g., molecules. The 
main intention of the work at hand is to investigate the 
transport behavior of finite quantum systems by means of 
an analysis based on the BE. To those ends we propose 
to identify the classical particle densities which appear 
in the BE with the the quantum mechanical occupation 
numbers of current eigenstates (|IV[) . To justify this con- 
cept we numerically and analytically analyze in Sect. [V] 
whether the dynamics of the current occupation numbers 
are indeed well described by an appropriate linear BE. 
This turns out to be the case but only for a quite specific 
form of the interactions within the model. For this type 
on interactions we concretely compute the adequate lin- 
ear BE and determine the diffusion coefficient using the 
form given in Sect. HH This coefficient turns out to be 
in accord with the results from the different approaches 
to similar systems mentioned above 0, [2l|, [22| • An ab 
initio numerical analysis of the full dynamics of the quan- 
tum system (including lattice site occupation numbers) 
shows that it indeed exhibits diffusive behavior controlled 
by the above diffusion coefficient. In the last Sect, we 
discuss the dependence of those results on special prop- 
erties of our model. 



II. BOLTZMANN EQUATION AND DIFFUSIVE 
SOLUTIONS 



As wellknown, in 1872 Boltzmann undertook to ex- 
plain the macroscopic dynamics of dilute gases. For 
the description of a gas Boltzmann introduced the \x- 
space, which is essentially a one-particle phase space. 
An A^-particle gas would thus technically be represented 
by N points in /i-space rather than one point in stan- 
dard Hamiltonian phase space. But instead of using TV 
points in /i-space for the description of the gas, Boltz- 
mann introduced a distribution function $(r,v,<) in a 
somewhat " coarse-grained" /i-space which is supposed to 
give the number of particles being in a -space cell around 
d 3 rd 3 v:=dxdydzdv x dv y dv z . Instead of trying to describe 
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the motion of every single particle (which is impossible 
due to the huge numbers of particles in a gas) Boltzmann 
suggested his famous equation which describes the time 
evolution of $(r, v, t) in /j,-space and is, in the absence of 
any external force, given by: 

— $(r, v) = v V r $(r, v) + $(r, v) scatt , (1) 

The first expression on the right-hand-side is supposed 
to account for the dynamics due to particles that do not 
collide, whereas <£>(r, v) scatt describes the dynamics aris- 
ing from collisions. Those dynamics are only taken into 
account in terms of the transition rates R and the (coarse 
grained) particle densities <&, neglecting all correlations 
and structures on a finer scale. Thus this type of dy- 
namics implement the famous "assumption of molecu- 
lar chaos" or the so-called "Stofizahlansatz" . (According 
to the Stofizahlansatz particles are not correlated before 
collisions, even though they get correlated by collisions, 
due to very many intermediate collisions before the same 
particles collide again). However, this treatment of col- 
lisions introduces the irreversibility into the BE that is 
not present in the underlying Hamiltonian equations. 
If either <E> is close to equilibrium or for systems in which 
particels only collide with external scattering centers (no 
particle-particle collisions), the BE takes on the following 
linear form. 



$(r, v) + v • V r $(r, v) 



#(v,v')$(r,v')dv', (2) 



A "velocity discretized" version of this linear BE reads 
*i(r) +v, • V r $,(r) =Y,Rij$j(r), (3) 

3 

where the matrix of rates Rij is of the standard form as 
appearing in master equations. (It is this discrete linear 
BE which we are eventually going to use to describe our 
quantum system.) 

As wellknown, in the limit of small density gradients 
("long wavelength hydrodynamic modes") ([3]) may fea- 
ture diffusive solutions, i.e., solutions that fulfill 

p(r) = «Ap(r). (4) 

where the particle density p(r) is given by 

p(r)=J2$i(*) (5) 

i 

The diffusion coefficient in ([4j) takes on the form 

K =-J2v i R^ 1 v j $° j . (6) 



$° is the equilibrium velocity distribution (i.e., the so- 
lution of Rij$® = 0, accordingly called "null-space"). 
RJa denotes the inverse of the rate matrix but without 
the null-space 0, [H| . (In the absence of any external 
scatterers k may diverge, indicating that no diffusive ba- 
havior can be expected.) 
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FIG. 1: Sketch of the systems which are investigated for trans- 
port: N identical weakly coupled subsystems with a non- de- 
generate ground state and a band of equidistant energy levels. 



III. DEFINITION OF THE MODEL AND ITS 
DIFFUSIVE BEHAVIOR 

The systems we investigate in the following are called 
"finite modular quantum systems" . (This type of system 
has also been investigated in[2(| [U, HH). Those sys- 
tems are essentially meant to allow for an investigation 
of transport behavior from first principles rather than 
to model some material or real physical system in great 
detail. Their total Hamiltonians may be written as 



JV 



N 



H = ^2h((j,) + ^2V((j,,p,+ l), with N + l=l 

(7) 

where h(p) denotes the local Hamiltonian of a subunit 
p, V(p, p + 1) a next-neighbor interaction between sub- 
units of the ring and N the total number of subunits in 
the system. A model which can be described by such 
a Hamiltonian is illustrated in Fig. Q] It consists of N 
identical subsystems where each subunit features a non- 
degenerate ground state, a wide energy gap (AE) and 
an energy band (Se) which contains n equidistant energy 
levels. The local Hamiltonian is defined by 



n+l 



5e 



Mm) = ^ EiaUp)ai(p) with Et=E + AE + i—, 

' J r? 



for i ^ 



(8) 



and <Ji(p) — \0,p)(i 1 p\. Here \0(i),p) denotes the ground 
state (i'th level of the excitation band) of the p'th sub- 
unit. The interaction is specified by 

n 

V(p,p + 1) = A %<5j(/i) ® d-jit* + 1) + h.c. (9) 

(h.c. denotes the hermitian conjugate of the previous 
sum). The Vij are complex numbers which are assumed 
to be adequately normalized such that A controlls the 
overall interaction strength. The choice of these Vij con- 
cretely determines the model. In [2^, [2l|, the Vij are 
simply chosen to be independent random numbers (with 
mean zero) which leads to diffusive dynamics (in a sense 
explained below). In the paper at hand it turns out that 
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this choice does not yield dynamics which are in accord 
with a BE. In fact, we will explain in Sect. [V]that it is 
only a rather specific choice which produces this accor- 
dance. However, at this point we do not specify the Vij 
any further. 

Irrespective of the choice for the the system may 
be viewed as a very simplified model for, e.g., a chain 
of coupled molecules or quantum dots, etc. In this case 
the hopping of the excitation from one subunit to an- 
other corresponds to energy transport. It may as well 
be viewed as a tight-binding model for particles on a lat- 
tice (in second quantization). There are, however, many 
(n) bands (orbitals per site) but no particle-particle in- 
teraction in the sense of the Hubbard-model (cf. (J7J), 
©)■ Thus one may characterize this model as a "single- 
particle multi-band quantum wire" with possibly random 
(interband) hoppings. Nevertheless, due to the indepen- 
dence of , Ei of p these are systems without disorder 
in the sense of, say, Anderson [26| ■ (Note that the "total 

particle number" N := YlTp °\ (M^Cm) is conserved) 
We define diffusive transport for this model based on the 
evolution of expectation values q^t) — (^f (i)\q(fj,)\^! (t)) 
of local fractions of some (globally conserved) quantity 
like energy or particle number (|$(t)) denotes the sys- 
tems full state). I. e., one may consider energy trans- 
port: q(p) = h(p) or particle transport: q(p) = n(p) = 
^™ =1 d-](^)(Ti(/i)., etc. We will call the transport diffu- 
sive if those local expectation values obey 



q^{t) = k ■ (q^+i(t) + q^-i{t) - 2q fl (t)) , 



(10) 



which is a discrete version of the diffusion equation (U) . 
In this case obviously only the diffusion coefficient k re- 
mains to be computed. If "diffusivity" is supposed to be 
a property of a system (TIT)]) has to apply of course to the 
largest part of all initial states under consideration. 



IV. CHOICE OF QUASIPARTICLES FOR 
MODELLING TRANSPORT 



It is not a priory clear which quantum observable 
could play the role of a classical particle density in phase 
space when one attempts to map the quantum dynamics 
of some system onto a BE. A minimum requirement is 
surely the clear correspondence to some velocity 
In the context of heat transport through electrically iso- 
lating crystals one identifies the particle (phonon) density 
with the average occupation number of a phonon mode. 
The velocity is then extracted from the phonon disper- 
sion relation [T3|. However, the phonons are eigen- 
modes of a harmonic chain (spring- and-ball- model), and 
the scattering arises from a (weak) anharmonic part of 
the interaction. In our model the interactions do not 
allow for a decomposition into an harmonic and an an- 
harmonic part. Thus this scheme cannot be applied here. 
In the context of particle transport through systems of in- 
teracting particles in periodic lattices, the (quasi)particle 



density is basically identified with the occupation number 
of eigenmodes (Bloch waves) of the interaction-free parti- 
cle, e.g., the crystal electron. The velocity (group veloc- 
ity) is then computed from the corresponding dispersion 
relation . Since our model (from a "Hubbard" point of 
view) features no particle-particle interaction such a de- 
scription would result in a BE without any scattering at 
all (Rij = 0), hence in ballistic transport. This however 
is not in accord with our findings (see below) . 
Thus, here we suggest to identify the particle density 
with the occupation number of eigenmodes of a suitable 
current operator. To define a current operator on the 
basis of the transported quantity we consider the time 
evolution of the corresponding local operator q{p) at site 
\x which is given by the Heisenberg equation of motion 

0.&§Ezh 



dt 



dt 



H,q(n) 



(11) 



since the operators q(n) mentioned above are explicitly 
time independent. After inserting equation (JT]) and ap- 



plying the explicit form of q{p)(= 
d 



dt 



- 1,m),<7(m) 



h(p), h(fx)) we obtain 
V{p,H + l),<z(/x) 



(12) 



If conserved quantities q(fJ-) are considered, currents are 
routinely defined on the basis of the temporal change of 
the respective densities by means of a (discrete) continu- 
ity equation which reads for q(fi): 

^(m) = J'(m>M + 1) -i0« - 1,M) = -divj. (13) 

Comparing Eq. (fl"2]) with Eq. (fl"3]) this suggests the defi- 
nition of a local current operator 



j Q (^,/i + l) = - [V(p,^ + l),q{fi) 
whereas the total currentoperator J® is given by 



(14) 



7V-1 



JV-1 



J°=£j°Gu,M+l) = ^£ [F(/,,M + l),g(M)](15) 

/I— 1 ^1—1 

E.g., in the case of energy transport the current operator 
is 

N—l n 

J H = ^ ■ A E i^M ® Oj-O* + 1) + h -c- (16) 

Since the current is a product of velocity and density we 
assign a velocity to our quasiparticles by the relation 

N 



Jn 



with Q = ^2q(p), 



(17) 



where is an eigenvalue of the corresponding current 
operator. In this way the velocities that eventually ap- 
pear in ([5]) may be defined. 
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V. ANALYSIS OF THE DYNAMICS OF THE 
MODEL 

In this Sect., we analyze whether the above described 
choice of quasiparticles is in accord with a BE from a 
dynamical point of view. Or, to formulate concise ques- 
tions: May the dynamics of the populations of the current 
eigenmodes as resulting from the Schrodinger-dynamics 
of the quantum model be described in terms of an ad- 
equate BE? And if so what would be the rates in the 
scattering term? 

In Sect. |lV]we suggested to identify the particle density 
as appearing in the BE by a quantity that corresponds to 
a specific velocity but not to a spatial coordinate. This 
is obviously in some sense insufficient since particle den- 
sities in phase space are labeled by velocity and poisi- 
tion. However, since the model features translational 
invariance, the current eigenmodes (the populations of 
which are supposed to correspond to the particle den- 
sities) stretch uniformly over the full model. Thus the 
dynamics of their populations may be expected to pos- 
sibly correspond to the dynamic as resulting from a BE 
for particle densities that are uniform with respect to the 
position coordinate, i.e., V r <I> = 0. In this case particle 
densities are only labled by velocities and may directly 
be identified with current eigenstate populations. In the 
above mentioned case, © simplifies to 

^n(t) =J2 R nm^>m(t). (18) 

m 

The above equation (fT5)) yields exponential decay for the 
$'s (possibly with various relaxation times). In the fol- 
lowing we investigate whether the same behavior results 
from the SE for the current eigenmode populations i.e., if 

we compute 3? n (t) from the definition $ n :=Tr//3(i)P n j 

where P n — \j n )(jn\ is the projector onto the subspace 
spanned by the current eigenstate \ j n ). If this is the case 
the SE and the BE may be in accord. However, this 
question is hard to answer in general without using nu- 
merical reasoning, since otherwise the quantum dynamics 
for 3> n (t) cannot be found. Thus, rather than analyz- 
ing the dynamics of the <&„(£) 's themselves, we analyze 
a function of those that can be estimated without us- 
ing numerics. Strictly speaking this of course means we 
go from a proof to a check of consistency (However, we 
also check dynamics of the $„ also directly numerically) . 
This function we call C(t) and construct it as 

C(t) =Y,3n ■ Tr [p{t)P n ] = Tr [j p(t)} , (19) 

n 

It is obviously just a weighted sum of the $„(£)'s. If 
its dynamics are in accord with (fT8)l . C(t) should also 
decay exponentially. As (fill))) shows, C(t) is simply the 
current expect at ionvalue. In the Heisenberg picture the 

latter reads C{t) =Tr| J(t)p(0)\. Again, we cannot an- 
alyze this in full generality, thus we specialize to a con- 
crete initial state, which is sometimes called a "deviation 



density matrix". It is given by p(0) = d^ 1 ■ 1 + eJ(0) 
(1 = identity, d = dimension of the corresponding (sub- 
space) and fulfills the relation Tr{/5o} = 1 for density 
operators due to the fact that the current operator is 

traceless. Thus we obtain C(t) = e-Tr{«7(i)J(0)}, which 

is simply the current-autocorrelation function. Without 
going into any detail here we should mention that, follow- 
ing concepts based on the Hilbert space average method 
(HAM) as presented in [l9|, |2l|, |23| , C(t) can be expected 
to reasonably describe the evolution of the expectation 
value of the current for almost any initial state. Thus 
the results which will be derived analytically below can 
safely be expected to apply to a much larger class of ini- 
tial states than covered by the deviation density matrix. 
Especially the results can be expected to apply to the 
largest part of all pure states, which is the class of states 
which will be primarily analyzed numerically below. 
The current-autocorrelation function reads: 

n 

C(t)= |(a|i(0)|/3}| 2 e-s( E *--^)-* ; (20) 

a,/3=l 

where |a(/3)), -E Q (/3) are energy eigenvectors respectively 
eigenvalues of the full, coupled system. 
Here and in the following we restrict ourselves to the "one 
excitation" (one-particle) subspace. This is possible since 
the particle number is conserved, cf. ([9]). If the coupling 
(A) is weak, it may be reasonable to approximate the 
true eigenvectors/eigenvalues of full system |a(/3)), E a (p) 
that appear explicitly in the correlation function, by the 
eigenvectors/eigenvalues of the uncoupled system from 
the one-particle subspace which feature the particle at a 
given site p. Since the current operator only "couples" 
states featuring the particle in adjacent sites, the double 
sum over sites collapses and we find in this approximation 

JV n 

C(t) « E E \(i,p\J(0)\»+lJ)\ 2 e-^ Bi - E ^. 

(21) 

For J(0) we plug in the energycurrent operator as given 
by (|T6|) . Obviously the addends do not depend on fi, thus 
performing the corresponding sum simply results in a 
prefactor N. If we assume AE >> Se and thus Ei » AE 
for the current operator (not for the exponential)we get 

n 

c~(*)=7E ifyiv-*^-*)-* 

(22) 

where 7 := 2N (-|Ai?) 2 (For particle transport we simply 
have to set AE = 1). In order to evaluate this expression 
we split up the double sum into two double sums. In the 
first double sum we perform the index transformation 
i = k — I + l,j = n — I + 1, in the second the index 
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transformation i = n — I + l,j = k — I + X. Thus using 
for Vij the definition from (f2"5)) yields: 



C{t) = 7 J ^A(fc,n)e~^^ (fe -" ) - t 
lfe=i 

n-1 ~j 

+ ^A'(fc,n)ei^( fc -")-* L , (23) 



with A(fc,n) = |Vfc-; + i,n-; + i| 2 

fe 

and A'(k,n) = ^ |K_ m , fe _; +1 | 2 . (24) 



Hence C(i) is essentially the Fourier transform of the A, 
A' . As explained above, if the quantum dynamics are 
claimed to be in accord with the BE, C(t) must decay 
exponentially. But this will only be the case if A, A' take 
the form of some Lorentzian in the argument (k — n). If, 
however, the Vy are chosen to be independent (gaussian) 
random numbers as done in [20|, the A, A' will simply be 
proportional to k and thus no exponential decay of C(t) 
is predicted within the framework of this approach. This 
expectation is confirmed by the numerical computation of 
C(t) as resulting from the Schrodinger equation (see Fig.'s 
[2J [3]) . Thus, in general, the dynamics of the quasiparticles 
(at least for this definition of quasiparticles) cannot be 
claimed to be in accord with a BE. If one modifies the 
weights of the Vij , however, one can enforce a Lorentzian 



shape upon A, A' . We now choose Vij as 



Va 



V, 



7(l-^-j|)-(l + « 2 (*-j) 2 )' 



(25) 



where Vij are still randomly distributed complex numbers 
normalized to n~ 2 J2ij \ Vij\ 2 = 1 an d a is arL (to some 
extend) arbitrary parameter. This choice yields in good 
approximation for large n: A w A' w n-(l — a 2 (k~ n) 2 ) -1 
and thus leads to 

7 ' I ^ 1 + a 2 (k ~n) 2+ ^[l + a 2 (k - nf 



(26) 



If a << 1 we may replace the sums by the following 
integrals 



C(t) 

TV) 



e i 



1 + a 2 (k — ny 



T dk 



1 + a 2 (k — nY 



rdk 




100 150 200 250 300 

time t (h/AE) 



400 



FIG. 2: Expectationvalue of the currentoperator obtained 
from the exact solution of the SE (here |i/>(0)) is a random 
pure state for which only the 2625 current eigenstates with 
the largest eigenvalues feature non-zero amplitudes) (dashed 
line). This is compared with the exponential function (solid 
line) featuring the relaxation time tr one obtains from an an- 
alytical analysis of C(t) (cf. (|20|) h System parameters: N = 
3, n = 3500, AE = 1.0, Se = 0.35, A = 5 ■ 10" 5 , a = 5 • 10" 3 . 
The pointed curve shows the current expectationvalue of a 
model without the proposed weighting function (cf. ([25}) the 
other system parameters are the same as mentioned above. 
The current relaxes, but not in an exponential way (see also 
Fig. 0. 



If furthermore (an) 2 >> 1 we may take the lower limit 
to negative infinity, obtaining 



C(t) w 2n 7 



cos[^(fc-rcH 
1 + o?(k — n) 



dk = C(0)e~^ (28) 



with 



C(0) 



2ir\ 2 AE 2 nN 



nha 



(29) 



This is obviously an exponential decay and thus indicates 
the applicability of an adequate BE. To countercheck this 
result, i.e., the validity of the above approximations, we 
compute the time evolution of some current eigenstate 
occupation number by solving the time dependent SE. 
The result is shown in the following figure (Fig. [2]). In 
Fig. [2] we use for the theoretical curve the results from 
the above analysis of C(t). Obviously there is rather good 
agreement between theory and numerics. Due to the fact 
that our theory just predicts only one relaxation time 
tr we conclude that for our model the "relaxation time 
approximation" seems to be valid. This finding enables 
as to specify the appropriate matrix of scattering rates 
for our model: 



-( - 

TR 



1 



N 



(30) 



where Sij denotes Kronecker's delta. After all this anal- 
(27) ysis it is justified to state that the microscopic dynamics 
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time t (h/AE) time t {h/AE) 



FIG. 3: Logarithmic plot of the current expectationvahies 
shown in Fig. [2] The straight dashed line corresponds to the 
numerical result for the modified interaction, the parabolic 
dashed line to the numerical result for a purely random in- 
teraction (the solid line corresponds to theory). The purely 
random interaction does not yield any exponential decay of 
the current. 



of the current eigenmode populations are consistent with 
a BE-description as given by (fT5| with the above matrix 
of scattering rates Rij . The equilibrium state <I>° is spec- 
ified by J2j 'I'': = 0. Thus one finds $° = 1/nN. All 
other states which are "orthogonal" to the equilibrium 
state (%2j $j = 0) correspond to eigenvectors of Rij with 
the eigenvalue — 1/tr. Hence the inverse Rij without 
the null-space is simply given by = — tr. With 

those results we may eventually evaluate the transport 
coefficient k according to (4). Consequently we insert in 
(4) for v n the velocity of the quasiparticles cf. (fT"7|) for 
which are approximatively v n rs Here we exploit 

that (j n \H\j n ) « AE. Plugging now all the results into 
(4) yields: 



n-N ■ AE 2 

TR 



II 



■ N ■ AE 2 



E 

Tr 



(31) 



where we exploited the invariance of traces with respect 

to unitary transformations. Since Tr|j 2 | is identical 

with the current-autocorrelation function at time t = 0, 
i.e., C(0), we may use (|2T)1) to find 



(32) 



This is the diffusion coefficient one obtains through an 
analysis based on an BE. If it is inserted into (flC)]) a sim- 
ple equation of motion for the populations of the subunits 
results. To check whether the dynamics produced by (10) 
coincide with the dynamics of the populations of subunits 



FIG. 4: Evolution of the occupation number of the subunits 
(n M (i)) for a pure initial state featuring ni(0) = 1 (112,3(0) = 
0). Dashed lines correspond to numerical solutions of the SE, 
solid lines to theory, cf. (I10|l . (132[) . The system parameters 
are the same like in Fig. [2] 



(n M (i)) obtained from direct numerical integration of the 
SE, we computed both. The result is displayed in Fig. 
[4] Obviously the agreement is rather good. Thus two 
conclusions can be drawn: i.) The model indeed shows 
diffusive behavior ii.) The diffusive behavior may be in- 
terpreted in terms of scattering quasiparticles. The scat- 
tering has to be treated as proceeding in such a way that 
the assumption of molecular chaos or the Stofizahlansatz 
apply. As a consequence, the macroscopic dynamics may 
be computed from a BE. 



VI. SUMMARY AND DISCUSSION 

We mainly demonstrated that the dynamics of a spe- 
cial class of finite modular quantum systems are to some 
extend in accord with the dynamics generated by an ad- 
equately set up BE. The BE here essentially appears as 
a rate equation rather than as an evolution equation for 
the phase-space density. The occupation numbers of cur- 
rent eigenmodes of the quantum system could be shown 
to obey this rate equation. Furthermore the occupation 
numbers of the local subunits ( "modules" ) of the quan- 
tum system evolve diffusively, with exactly the same dif- 
fusion coefficient that one gets from analyzing the be- 
havior of long-wavelength hydrodynamical modes of the 
BE. 

However, we consider it crucial, that the above de- 
scribed applicability of a BE to modular quantum sys- 
tems does not hold in general. In the example consid- 
ered in this paper the applicability has been enforced by 
the special form of the interaction as described in (|2"5j) . 
This special form contains no restriction regarding the 
phases of the transition (interaction) matrix elements 
but requires that there weights essentially fall off in a 
Lorentzian shape with energy differences getting larger. 
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A statement about a the "typicality" of such interactions 
can hardly be made, but as a mathematical condition the 
special form appears quite restrictive. In contrary to this 
large classes of those modular quantum systems exhibit 
diffusive behavior with respect to occupation numbers 
of their subunits even if the interaction does not fea- 
ture the above form. This implys that there may be a 
(large?) class of systems in general that exhibit diffusive 
behavior but cannot be described in terms of a BE pic- 
ture, i.e., the concept of scattering quasiparticles may, 
strictly speaking, be inapplicable. Most investigations of 
transport based on the Kubo formula focus on the ques- 
tion whether the integral over the current-autocorrelation 
function is finite, not on whether the correlation function 
decays exponentially. But the latter would be needed for 
the applicability of a BE. 

Also the fact that a translationally invariant "one- 
particle" system exhibits diffusive behavior requires ex- 
planation. According to standard solid-state theory exci- 
tations should correspond to free "lattice-particles" fea- 
turing a dispersion relation depending on the periodic 
Hamiltonian. Thus transport is expected to be ballis- 
tic. However, in the case of a rather small amount of 
subunits and a large amount of "orbitals" per subunit 
the band structure looks more like a disconnected set 
of points in the E(k) vs. k diagram, rather than the 



usual set of smooth dispersion relations. In the same case 
eigenstates of the current operator do not coincide with 
Bloch energy eigenstates. They only do coincide in the 
limit of the number of subunits going to infinity in which 
one then also gets smooth dispersion relations. So for a 
comparatively small number of subsystems the current 
eigenstates are no stationary states. Thus it appears to 
be especially the limit of small numbers of subunits that 
yields the diffusive behavior. (We expect that to vanish 
in the limit of infinitely many subunits, investigations in 
that direction are currently being done.) Thus, for the 
case of a small number of subunits we have the seem- 
ingly paradox situation that a translationally invariant 
one particle system maps onto a BE in which diffusiv- 
ity only arises from external scatterers. This, however, 
one routinely only expects for systems featuring disorder, 
defects or impurities. 



Acknowledgments 

We thank K. Barwinkel and J. Schnack for fruitful dis- 
cussions. Financial support by the Deutsche Forschungs- 
gemeinschaft and the Graduate College 695 "Nonlincar- 
ities of optical Materials" is greatfully acknowledged. 



R. Kubo, M. Toda, N. Hashitsume, Statistical Physics II. [15] 
Nonequilibrium Statistical Mechanics, Springer- Verlag, [16] 
Berlin (1991) [17] 
G. D. Mahan, Many Particle Systems, Plenum, New York 
(1990) 

A. Klumper and K. Sakai: J. Phys. A 35, 2173 (2002) [18] 
X. Zotos, F. Naef, P. Prelovsek, Phys. Rev. B. 55, 11029 
(1997) [19] 
F. Heidrich-Meisner, A. Honecker, D. Cabra, W. Brenig, [20] 
Phys. Rev. B. 66, 140406 (2002) 

P. Jung, R. W. Helmes, A. Rosch, Phys. Rev. Lett. 96, [21] 
067202 (2006) 

L. Boltzmann, Lectures on Gas Theory, University of [22] 
California Press, Los Angeles (1964) 

C. Cercignani, The Boltzmann Equation and Its Applica- [23] 
tions, Springer, N. Y. (1988) 

R. E. Peierls, Quantum Theory of Solids, Claredon Press, [24] 
Oxford (2001) 

L. P. Kadanoff, G. Baym, Quantum Statistical Mechan- 
ics, Benjamin, New York (1962) [25] 
J. Rammer, H. Smith, Rev. Mod. Phys. 58, 323 (1968) 
K. Barwinkel, Z. Naturforsch. A 24A, 22 (1969); 24A, [26] 
38 (1969) [27] 
W. Kohn, J. M. Luttinger, Phys. Rev. 108, 590 (1957) 
K. Aoki, J. Lukkarinen, H. Spohn, J. Stat. Phys. 124 
(2006) 



K. Hornberger, Phys. Rev. Lett. 97, 060601 (2006) 
B. Vacchini, Int. J. Theor. Phys. 44, 1011 (2005) 
J. Gemmer, M. Michel, G. Mahler Quantum Thermody- 
namics, Lecture Notes in Physics 657, Springer- Verlag, 
Berlin (2004) 

H. P. Breuer, F. Petruccione, The Theory of Open Quan- 
tum Systems, Oxford University Press, Oxford (2002) 
J. Gemmer, M. Michel, Eur. Phys. J. B 53, 517 (2006) 
M. Michel, J. Gemmer, G. Mahler, Phys. Rev. Lett. 95, 
180602 (2005) 

J. Gemmer, R. Steinigeweg, M. Michel, Phys. Rev. B, 
73, 104302 (2006) 

R. Steinigeweg, J. Gemmer, M. Michel, Europhys. Lett. 
75, 406 (2006) 

H. P. Breuer, J. Gemmer, M. Michel, Phys. Rev. E 73, 
016139 (2006) 

R. Balescu, Equlibrium and Nonequilibrium Statistical 
Mechanics, John Wiley & Sons, New York, London, Syd- 
ney, Toronto (1975) 

W. Brenig, Statistical Theory of Heat, Springer- Verlag, 

Berlin, Heidelberg, New York (1989) 

P. W. Anderson, Phys. Rev. 109, 1492 (1958) 

M. Michel, J. Gemmer, G. Mahler, Euro. Phys. J. B 42, 

555 (2004) 



